ggplot2

library(ggplot2)
ggplot(diamonds[sample(nrow(diamonds), 1000), ], aes(carat, price)) +
  geom_point(alpha = 0.2) +
  geom_smooth() +
  facet_wrap(~cut) + ggtitle("Diamond price by carat and cut")

Lattice

lattice::show.settings()

Base

pairs(
  iris[,1:4], pch = 19, 
  col = iris$Species, lower.panel = NULL
)